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ABSTRACT 


A  Fortran  subroutine  to  generate  all  or  particular  combi¬ 
nations  of  n  things  taken  m  at  a  time  is  herein  described. 
The  procedure  establishes  and  tabulates  a  certain  one-to-one 
correspondence  between  the  integers  from  1  to  nCm  and  the 

individual  combinations  of  n  distinct  "things"  taken 
m  at  a  time,  where 


C  =  n! 
n  m  mTCn-my! 

In  all  cases,  the  set  of  "things"  treated  by  the  subroutine 
is  the  set  of  integers  {l,...,n}. 


PROBLEM  STATUS 


This  is  a  final  report  on  this  phase  of  the  problem. 


AUTHORIZATION 


NRL  General  and  Administrative  Function  78-1601 
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NRL  Memorandum  Report  2340 
NRL  Computer  Bulletin  25 


A  Fortran  Subroutine  to  Produce  Combinations 
of  n  Distinct  Things  Taken  m  at  a  Time 

1.0  IDENTIFICATION 

1.1  Title 

Combinations  Subroutine 

1 . 2  Identification  Name 
GO-NRL-CMBN 

1 . 3  Classification  Code 

GO  -  Statistical  Analysis  and  Probability,  General 

1 . 4  RCC  Identification  Name 
G0001R00 

1 . 5  Entry  Points 
CMBN 

1.6  Programming  Language 

Language:  3600/3800  F0RTRAN 

Routine  Type:  Subroutine 
Operating  System:  Drum  Scope  2.0 

1 . 7  Computer  and  Configuration 
CDC  3800 

1 . 8  Contributor  or  Programmer 

James  P.  Grimes,  Code  7816JG,  Research  Computation 
Center,  Mathematics  and  Information  Sciences 
Division 

1 . 9  Contributing  Organization 

NRL  -  Naval  Research  Laboratory,  Washington,  D.  C. 
20390 
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1.10  Program  Availability 

1.10.1  Submittal:  Program  write -up ,  Fortran  source 

deck,  source  listing 

1.10.2  On  File:  RCC  Program  Library 

1.11  Verification 


The  subroutine  was 

checked  for  the  following 

values  of  N,  M,  and 

NN; 

(see  Section  3.2  for 

definition  of  N,  M, 

and 

NN)  . 

N 

M 

NN 

17 

15 

0 

17 

16 

0 

20 

20 

0 

101 

100 

2 

0 

0 

1 

20 

20 

1 

19 

0 

0 

1 

1 

0 

1 

1 

1 

6 

7 

0 

0 

1 

0 

1.12  Date 

1  July  1971 
2.0  PURPOSE 

2.1  Description  of  the  Routine 

To  generate  all  or  particular  combinations  of  the 
integers  from  1  to  n  taken  a  given  number  at  a  time. 

This  procedure  establishes  and  tabulates  a  certain 

one-to-one  correspondence  between  the  integers 

from  1  to  C  and  the  individual  combinations  of 
n  m 

n  distinct  "things"  taken  m  at  a  time.  Here 


m ! (n-m) ! 

In  all  cases,  the  set  of  "things"  treated  by  the 
subroutine  is  the  set  of  integers  {l,...,n}. 
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The  type  of  correspondence  referred  to  may  be 
illustrated  for  the  case  n=5,  and  m=3,  as  follows: 


Index 


Combination 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


1.2.3 

1.2.4 

1.2.5 

1.3.4 

1.3.5 

1.4.5 

2.3.4 

2.3.5 

2.4.5 

3.4.5 


2 . 2  Problem  Background 

In  many  phases  of  investigative  mathematics, 
statistics,  probability  and  related  areas, 
the  need  often  arises  to  generate  and  tabulate 
all  or  particular  combinations  of  n  distinct 
"things"  taken  m  at  a  time.  A  computerized 
algorithm  is  the  only  practical  way  to  handle 
this  problem,  especially  as  the  number  of 
combinations  and  the  number  to  be  taken  at 
a  time  become  relatively  large.  Previously, 
this  type  of  computerized  procedure  was 
unavailable  to  NRL's  scientific  community. 


3.0  USAGE 

3 . 1  Calling  Sequence  or  Operational  Procedure 

The  calling  sequence  for  the  subroutine  is: 

CALL  CMBN  (N,  M,  NN,  MM) 

All  four  parameters  are  in  CDC  3600/3800  integer 
format. 

3 . 2  Arguments,  Parameters,  and/or  Initial  Conditions 
Input  parameters : 

N  and  M  are  the  n  and  m  described  in  Section  2.1, 

i . e . ,  as  in  „C  . 

’  n  m 

NN  is  an  input  integer  which  establishes  the 
correspondence  referred  to  in  Section  2.1. 

If  NN=0,  the  correspondence  is  set  by  the  subroutine 
and  all  combinations  from  1  to  C  are  generated, 
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indexed,  and  printed  out  by  the  subroutine.  Otherwise, 

if  NN  =  k,  k  =  1,...,  C  ,  only  that  combination 

n  m 

which  corresponds  to  k  is  generated,  indexed,  and 
printed  out  as  demonstrated  in  Section  2.1. 

MM  is  a  one-dimensional  input  array  which  must 
be  dimensioned  in  the  calling  program  and  whose 
length  must  be  at  least  4(m+l). 

Output  parameters: 

There  are  no  output  parameters  as  such  since  all 
output  is  performed  by  the  subroutine.  However, 
after  execution  of  the  subroutine  for  each  value 
of  NN,  whether  supplied  by  the  user  or  the  subroutine 
itself,  the  m  integers  comprising  the  output 
combination  occupy  the  first  m  locations  of  the 
array  MM.  The  subroutine  prints  out  this  combination 
and  the  MM  array  is  reused  for  the  next  combination. 
This  procedure  is  repeated  until  all  combinations 
(NN  =  0)  or  the  particular  combination 
(NN  =  k,  k  =  l,...,nC  )  is  generated  and  printed  out. 

3.2.1  Special  Provisions  and  Restrictions 

Apparently,  the  basic  algorithm  imposes 

no  restrictions  on  the  parameters.  However, 

there  are  computer-imposed  limitations, 

such  as  storage  capacity,  word  size,  etc. 

Nevertheless,  to  facilitate  the  output, 

n,  m,  and  C  have  been  limited  to 
^  n  in 

10s.  Of  course,  n>  m. 

For  those  cases  in  which  m  £  16,  CMBN  uses 

logarithms  for  some  of  the  intermediate 

steps;  this  fact  results  in  appreciably 

longer  computation  time  for  these  cases. 

For  the  most  part,  computation  time  depends 

on  the  size  of  m  and  C  . 

n  m 

3 . 3  Space  Required  (Decimal  and  Octal) . 

3.3.1  Unique  Storage:  1475  octal  (829  decimal) 
locations  m  CMBN 

The  dimension  of  the  input  array  MM  is  user 
supplied . 

3.3.2  Common  Blocks :  None 

3.3.3  Temporary  Storage :  None 
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3.4  Messages  and  Instructions  to  the  Operator 
None 

3 . 5  Error  Returns,  Messages,  and  Codes 
None 

3 . 6  Informative  Messages  to  the  User 

The  message  "INVALID  INPUT  PARAMETER (S)  N,M*0  OR 

M>N"will  be  printed,  if  the  input  parameters  are 
invalid. 


3 . 7  Input 

See  the  "Input  parameters"  described  in  Section  3.2 

3 . 8  Output 

See  "Output  parameters"  in  Section  3.2 

3 . 9  Formats 
None 

3.10  External  Routines  and  Symbols 
None 


3.11  Timing 


Running  time  is  primarily  dependent  on  the  size  of 

m  and  C  .  Some  examples  of  running  time  on  the 
n  m  r  6 


CDC-3800  computer 
24  seconds;  i5^5> 


follow:  7C4 ,  24  seconds;  gC^, 

72  seconds;  20^17’  ^  secon^s- 


3.12  Accuracy 


The  procedure  has  been  tested  by  the  author  for 
values  of  m  up  to  100.  Apparently,  it  continues 
to  be  valid  for  higher  values  of  m,  but  computation 
time  greatly  increases. 


3.13  Cautions  to  Users 

3.13.1  MM  is  an  input  array  and  must  be  dimensioned 
in  the  main  program.  Its  size  must  be  at 
least  4(m+l).  The  array  name  in  the  main 
program  may  be  any  integer  variable. 
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3.13.2  Since  the  value  of  C  can  become  very 

n  m 

large,  the  user  should  make  sure  he  does 
not  exceed  storage  capacity  or  machine 
word  size  using  integer  format.  To 
facilitate  the  output  by  the  subroutine 
n,  m,  and  C  must  be  less  than  10^. 

3.13.3  The  user  should  consider  running  time  and 

the  amount  of  output  produced  for  large 

values  of  m  and  C  .  • 

n  m 

3.13.4  Of  course,  n  and  m  must  be  positive 
and  n  >  m . 

3.14  Program  Deck  Structure 
gJ0B  card 

g FTN '  card 

main  program 
CMBN  subroutine 

SC0PE 
gL0AD  card 

gRUN  card 

gg (end-of -f ile  card) 

3.15  References  -  Literature  -  Appendices 

Narec  Bulletin  #120,  Neliac  Bulletin 
#58,  by  Peter  C.  Ryan,  formerly  Code 
6410,  NRL . 

4.0  METHOD  OR  ALGORITHM 

The  procedure  used  in  the  subroutine  is  a  modification 
and  generalization  of  an  algorithm  used  in  a  NELIAC-N 
computer  function  by  Peter  C.  Ryan,  formerly 
Code  6410,  NRL.  For  further  information  or  details, 
see  reference  listed  in  Section  3.15. 
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o-i 

O'*: 


nr, 

c~ 

5.0  SOURCE  LANGUAGE  LISTING 


SUBROUTINE  CMBN (N*M»NN»MM)  000 

00-NRL-CM9NiCDC  3800  COMPUTER  SUBROUTINE  TO  PRODUCE  100 

COMBINATIONS  of  n  DISTINCT  things  TAKEN  M  at  a  TIME,  200 

(IN  ALL  CASES*  THE  SET  OP  THINGS  TREATED  BY  THE  300 

SUBROUTINE  IS  THE  SET  OP  INTEGERS  *1 » *  » . *N=)  *00 

TITLE-COMBINATIONS  SUBROUTINE  500 

IDENT  NAME-GO-NRL-CMBN  600 

LANGUAGE-PORTRAN  TOO 

PROGPAMMER-JAMES  P,  GRIMES, CODE  7816JG  BOO 

RESEARCH  COMPUTATION  CENTER*  MlSo  900 

(A  MODIPICATION  AND  GENERALIZATION  Op  A  1000 

NELIAC-N  COMPUTER  punction  BY  MR,  noo 

PETER  C.  RYAN*PORMERLY  CODE  6410, nRL)  1200 

IDENT  NUMBER-GOOOlROO  1300 

ORGANIZATION. NRL-NAVAL  RESEARCH  LABORATORY  1400 

WASHINGTON»D,C, 20390  1500 

DATE-2  NOVEMBER  1970  1600 

PURPOSE-TO  GENERATE  ALL  OR  PARTICULAR  COMBINATIONS  1700 

op  the  integers  prom  i  to  n,  iboo 

1900 

U9AGE-THE  CALLING  SEQUENCE  FOR  THE  SUBROUTINE  IS  I  2000 

CALL  CMBN(N*M*NN,MM)  2100 

ALL  POUR  PARAMETERS  ARE  IN  CDC  3600/3800  2200 

INTEGER  FORMAT.  2300 

2400 

INPUT  PARAMETERS-N  AND  M  ARE  THE  INPUT  INTEGERS  IN  2500 

C(N,M),  THAT  IS,  THE  COMBINATIONS  2600 

OP  N  THINGS  TAKEN  M  AT  A  TIME*  2700 

2800 

NN, AN  INPUT  INTEGER  WHICH  ESTABLISHES  2900 

A  CERTAIN  ONE-tO-ONE  CORRESPONDENCE  3000 

BETWEEN  THE  INTEGERS  PROM  1  TO  C(N,M)  3100 

AND  THE  INDIVIDUAL  COMBINATIONS  OF  N  3200 

things  taken  M  at  A  TIME.  3300 

3400 

IF  NNaO,  THE  CORRESPONDENCE  IS  SET  BY  3500 

THE  SUBROUTINE  AND  ALL  COMBINATIONS  3600 

FROM  1  TO  C(N,M)  ARE  GENERATED,  3700 

INDEXED*  AND  PRINTED  OUT  BY  THE  3800 

SUBROUTINE,  3900 

4000 

IP  NNaK,  K«1,...,C(N»M) ,  ONL*  THAT  *100 

COMBINATION  WHICH  CORRESPONDS  TO  K  *200 

IS  GENERATED,  INDEXED,  AND  PRINTED  *300 

OUT.  **00 

4500 

MM  IS  A  ONE-DIMENSIONAL  INPUT  ARRAY  *600 

WHICH  MUST  BE  DIMENSIONED  IN  THE  *700 

CALLING  PROGRAM  AND  WHOSE  LENGTH  *800 

MUST  BE  AT  LEAST  4*(M*1>.  *900 
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output  parameters-there  are  no  output  parameters 

AS  SUCH  SINCE  ALL  OUTPUT  IS 
PERFORMED  BY  THE  SUBROUTINE. 

RESTRICT I ONS-TO  FACILITATE  THE  OUTPUT,  N.M.AND 
C (N.M)  HAVE  BEEN  LIMITED  TO  10#«5. 

N  AND  M  MUST  BE  POSITIVE  AND  N2M, 

COMPUTATION  TIME  DEPENDS.  FOR  THE  MOST  PART,  ON 
THE  SIZE  OF  M  AND  C<N,M>. 

THE  PROCEDURE  HAS  BEEN  TESTED  FOR  VALUES  OF  M  UP 
TO  100. 

CAUTION-AS  M  AND  C(N,M>  BECOME  LARGE,  THE  RUNNING 
TIME  AND  AMOUNT  OF  OUTPUT  INCREASES 

appreciably. 


DIMENSION  JSPEC (5) 

DIMENSION  MM ( 1 > 

IF<M,LE.0.OR.M.GT.N)60,70 
60  PRINT  65 

65  FORMAT (5X,*INVALID  INPUT  PARAMETER(S)  N.MjO  OR  M>N*,3(/)) 

return 

TO  IUKpO 

IF (NN.EQ.O) 1,14 

1  I*IK«i 
IF(N,EQ.1)2,3 

2  FNF«1, 

GO  TO  5 

3  FNF*1, 

DO  4 » 1*1 »N 

4  ENF»FNF*I 

5  IF (M.EQ, 1 ) 6,7 

6  FMF"1. 

GO  TO  9 

7  FMF*1 , 

DO  8» 1*1 »M 

8  FMF*FMF*I 

9  K*N-M 

IF<K. EQ.O. OR. K.EQ.l) 10,11 

10  FKF»1 • 

GO  TO  13 

11  FKF*1 . 

DO  12, 1*1, K 

12  FKF*FKF*I 

13  NC»FNF/(FMF*FKF>*.5 
DO  49.  11*1, NC 


5000 

5100 

5200 

5300 

5400 

5500 

5600 

5700 

5800 

5900 

6000 

6100 

6200 

6300 

6400 

6500 

6600 

6700 

6800 

6900 

7000 

7100 

7110 

7120 

7130 

7140 

7150 

7200 

7300 

7400 

7500 

7600 

7700 

7800 

7900 

8000 

8100 

8200 

8300 

8400 

8500 

8600 

8700 

8800 

8900 

9000 

9100 

9200 

9300 

9400 
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NN-II 

9500 

14 

MN"4* (M*l ) 

9600 

DO  15  e I"1 *MN 

9700 

15 

MMtI>aO 

9800 

J  MaM* } 

9900 

ie"IM*M*l 

10000 

IN»IC*M* 1 

10100 

IF (M.FQ.l) 16.17 

10200 

16 

MM»NN 

10300 

00  TO  45 

10400 

17 

MM ( I C ♦ 1 ) *N 

10500 

MM ( IN* 1 ) »M 

10600 

Il>0 

10700 

IF(M.6T. 15)30*18 

10800 

18 

MlaM-1 

10900 

DO  29» J»1 *Ml 

11000 

MM(IC*J*1)"MM(IC*J) »MM ( IN* J) *1 

moo 

1 a*mm ( IC* J* 1 ) 

11200 

IBal 

11300 

IH«IL 

11400 

MJaM-J 

11500 

DO  19,L«1»MJ 

11600 

19 

TB*IB*l 

11700 

DO  26»I»1»IA 

11800 

IO"l 

11900 

DO  20 • L“1 »M J 

12000 

20 

TD-ID*(IA-I*L> 

12100 

IS-ID/IB 

12200 

Il-IH 

12300 

IH-Il*IE 

12400 

IF(IH,GT,NN)21#22 

12500 

21 

MM(IN*J*l)aI 

12600 

GO  TO  27 

12700 

22 

IF(NN,EQ,IH)23i26 

12800 

23 

MM ( IN* J*1 ) ■! 

12900 

IF«0 

13000 

DO  24*L“1 * J 

13100 

24 

IF»IF*MM ( IN*L*1 > 

13200 

MM(IN*J*2)«N-M-IF*J*1 

13300 

J2a J*2 

13400 

DO  25»L"J2*M 

13500 

25 

MM(IN*|_*1)»1 

13600 

GO  TO  42 

13700 

26 

CONTINUE 

13800 

27 

IF(M-1.E0.J>28»29 

13900 

28 

MM<IN*M*1)«NN-IL 

14000 

29 

continue 

14100 

GO  TO  42 

14200 

30 

MlaM-1 

14300 

DO  41  * J«1 »M1 

14400 
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14500 

IA*MM(IC*J*1) 

14600 

FB-0 

14700 

IH-IL 

14800 

M J«M- J 

14900 

00  31 *L*1 »MJ 

15000 

Fl»L 

15100 

FL«L0GF(FL> 

15200 

31 

F8«F6*FL 

15300 

DO  38* 1*1 t I A 

15400 

F0«0 

15500 

DO  32*L«1iMJ 

15600 

ID*IA-I*L 

15700 

FA»ID 

15800 

FA«LOSF (FA) 

15900 

32 

fo»fo*fa 

16000 

FOuFD-FB 

16100 

FO»EXP(FD) *.5 

16200 

IE"FO 

16300 

IL-IH 

16400 

IH»IL*IE 

16500 

IF ( IH.GT.NN) 33(34 

16600 

33 

MM(!N*J*1)*I 

16700 

GO  TO  39 

16800 

34 

IF(NN,EQ.IH)35*38 

16900 

35 

17000 

IF-0 

17100 

DO  36»L»1»J 

17200 

36 

TF«IF*MM(IN*L*1) 

17300 

MM(IN*J*2>«N-M«IF*J*1 

17400 

J2«J*2 

17500 

DO  37»L*J2*M 

17600 

37 

MM(IN*L*1)«1 

17700 

GO  TO  42 

17800 

36 

CONTINUE 

17900 

39 

IF  (M-l.EQ.  JU0.41 

18000 

40 

MM<IN*M*1)«NN-IL 

18100 

41 

CONTINUE 

18200 

42 

DO  43tK«l.M 

18300 

43 

MM(  IM*K*1 ) «HM (IM*K)*MM(IN*K*1) 

18400 

00  44*  XX*1 *M 

18500 

44 

MM(IX)«MM(IM*IX*1) 

18600 

IF (IJK»EG# 1)45*46 

18700 

45 

IF(NC.LT,10E4)  LL»5 

18800 

IF<Nc.LT.10E3)  U«6 

18900 

IF(NC.LT.10E2)  LL"3 

19000 

IF(NC.LT.lOEl)  UL»2 

19100 

IF (NC*LT . 10E0)  llml 

19200 

GO  TO  47 

19300 

46 

IF(NN.LT.10E4)  LL«5 

19400 

10 


47  IF (N»LT • 10E4)  JJ»5 
IF ( N,LT# 1 0E3 )  JJ*4 
IF (N.LT, 10F2)  JJ»3 
IF(N,tT.10Fl)  JJ»2 
IF(N.LT.IOEO)  JJ»1 
JKL»5*LL 

KK»(130-LL-5)/(JJ*D 

ENCODE (37*48* JSPEC) LL*KK*JJ*JKL*KK»JJ 

48  FORMAT(5H(3X*I*IIt4H,2x*l2,2H(I*lIt6H»lX)/<» 
C  I2*2HX*.I2»2H<I*I1*6H,1X> > > > 

PRINT  JSPEC*NN* (MM( J) • J«1*M) 

49  CONTINUE 
END 


19500 

19800 

19700 

19800 

19900 

20000 

20100 

20200 

20300 

20400 

20500 

20600 

20700 


11 


UNCLASSIFIED 


6.0 


COMPARISON 


No  routine  was  available  to  the  author  for  comparison. 
7.0  TEST  METHOD  AND  RESULTS 


A  test  program  that  demonstrates  all  combinations  of 
7  "things"  taken  4  at  a  time  and  specific  combinations 
in  accordance  with  the  established  correspondence  of 
the  individual  combinations  and  the  positive  integers 
along  with  the  subroutine  printouts  follow. 


PROGRAM  TCMBN 

DIMENSION  II (20) *  KARRAY (5) 
DATA  (KARRAY«3fX2»7»20»5) 
N«7 

M«4 

CALL  CMBN(N*M*0f II) 

PRINT  1 

1  FORMAT ( 10 (/) ) 

00  2  1*1 *5 

CALL  CMBN (N#M»KARRAY (I) *11) 

2  CONTINUE 
END 
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(Page  14  is  Blank) 

1  1  2  3  4 

2  12  3  5 

3  1  2  3  <S 

4  1  ?  3  7 

5  12  4  5 

4  12  4  4 

7  12  4  7 

3  12^4 

Q  12  4  7 

10  1247 

11  1344 

1?  1  .3  4  * 

13  13  4  7 

14  13  4  4 

14  13  4  7 

14  13  4  7 

17  1  4  4  f, 

14  14  4  7 

IQ  14  4  7 

20  1547 

21  2345 

22  2344 

23  2347 

24  2344 

25  2357 

24  2347 

27  2  4  4  4 

25  2457 

29  2447 

30  2  5  4  7 

31  3444 

32  3  4  5  7 

33  3447 

34  3547 

35  4  5  4  7 


3  12  3  4 

12  13  4  4 

7  12  4  7 

20  1547 

5  12  4  5 
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